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We perform a coupled-channel study of the possible deuteron-like molecules with two heavy flavor quarks, in- 
cluding the systems of £) ( *'D ( *> with double charm, with double bottom and D^'B^ with both charm and 
bottom, within the one-boson-exchange model. In our study, we take into account the S-D mixing which plays 
an important role in the formation of the loosely bound deuteron, and particularly, the coupled-channel effect in 
the flavor space. According to our calculation, the states D l " ) D i * ) [I(J p ) = 0(1 + )] and (D ( * ) D ( * ) ) s [7 /> = 1 + ] 
with double charm, the states B M B M [I(J P ) = 0(1 + ), 0(2 + ), 1(0 + ), 1(1 + ), 1(2 + )], (S w S w ) s [y /> = + ,l + ,2 + ] 
and (B M B M ) SS [J P = + ,l + ,2 + ] with double bottom, and the states D M B M [I(J P ) = 0(0 + ),0(l + )] and 
(D ( * ) B ( * , )j[y p = + , 1 + ] with both charm and bottom are good molecule candidates. However, the existence 
of the states D (,) D (t) [I(J p ) = 0(2+)] with double charm and D (,) B ( ' ) [I(J P ) = 1(1+)] with both charm and bottom 
is ruled out. 

PACS numbers: 14.40.Rt, 14.40.Lb, 12.39.Hg, 12.39.Pn 



I. INTRODUCTION 

More and more experimental observations have stimulated 
the extensive discussions of exotic states. The molecular state 
explanation to the reported charmonium-like states X, Y, Z be- 
comes popular due to the fact that many charmonium-like 
states near the threshold of charmed meson pair, i.e., 

X(3872) ~ m D n>, K(3930) ~ m D , D , 
7(4140) ~ m mBl , F(4274)~m Djo(2317)4 . 

And, in the past decade there is abundant literature with the 
study of the heavy flavor molecular states lflU23ll. 

The concept of molecular state with hidden charm was first 
proposed by Voloshin and Okun thirty years ago and they 
studied the interaction between the charmed and anti-charmed 
mesons JU Later, De Rujula, Georgi and Glashow sug- 
gested that the observed t/r(4040) is a D*D* molecule 12511 . 
By the quark-pion interaction model, Tornqvist investigated 
the possible deuteron-like two meson bound states with BB* 
or B*B* component I26ll27]l . At present, carrying out the phe- 
nomenological study of the heavy flavor molecular state is still 
a hot research topic of hadron physics. 

Usually, the hadron configurations mainly include 



Meson : qq, Qq, QQ 

Baryon : qqq, Qqq, QQq, 
Molecular state 



Hadron 



Exotic state : 



Hybrid 
Glueball 



where q and Q denote the light (u, d, s) and heavy (c, b) 
quarks, respectively. Among the conventional baryon states, 
the baryons with double charm or double bottom is of the QQq 
configuration. The SELEX Collaboration reported the first 
observation of a doubly charmed baryon E* c in its charged 
decay mode S* c — > A+A"~;r + fl28ll and confirmed it in the de- 
cay mode E* c — > pD + K~ 12911 . However, later the BABAR 
Collaboration searched for E+ c in the final states A^K~n + and 
Bj?7r + , and HJ C + in the final states A^K~n + n~ and Sj.V + 7t + , and 
found no evidence for the production of the doubly charmed 
baryons l30ll . The Belle Collaboration reported no evidence 
for the doubly charmed baryons in the final state A^K^n + , 
either [31]. Although these doubly charmed baryons were 
not confirmed by BABAR and BELLE, it is still an interest- 
ing research topic to search for such doubly charmed baryons 
experimentally. 

Besides the doubly heavy flavor baryons, it is also very in- 
teresting to study other systems with two heavy flavor quarks. 
The heavy flavor molecular state with two charm quarks pro- 
vides another approach to investigate the hadron states with 
double charm. For this kind of hadron, its typical configu- 
ration is [cg][«/]. To answer whether there exist such heavy 
flavor molecular states with double charm or not, in this paper 
we apply the one-boson-exchange (OBE) model to perform a 
dynamic calculation of their mass spectroscopy. This study 
is not only a natural extension of the previous work of the 
heavy flavor molecular state with hidden charm, but also pro- 
vides new insight into exploring the hadron states with dou- 
ble charm. Besides the hadron states with double charm, we 
also investigate the hadron states with double bottom and the 
hadron states with both charm and bottom. 
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This paper is organized as follows. After the introduction, 
we present the derivation of the effective potential in Sec- 
tion|IIl We summarize our numerical results and perform some 
analysis in Section [III] and draw some conclusions in Sec- 
tion[IV] We also give some useful formulas in the Appendix. 
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II. FORMALISM 



A. The lagrangians and the coupling constants 

In the present paper, we investigate the possible molecules 
of (Z)WZ)W) with double charm, (5WjjOO) with double bottom 
and (D^'B^*') with both charm and bottom. In our study, we 
take into account the S-D mixing which plays an important 
role in the formation of the loosely bound deuteron and, par- 
ticularly, the coupled-channel effects in the flavor space. We 
study the systems with total angular momentum J < 2. We 
list the channels for different systems in Tables H1HI1 

The Lagrangians under the heavy quark symmetry and the 
SU(3)-flavor symmetry read I32U35I1 



-Chhv 



HHcr - 



gs Tr[//(fiW?«»], 



where H iQ) and H (Q> are defined as 
1+ 



H (Q) 



and 



H¥>= 7 oH^ 7o = [P?y,+P a75 ] 



/(G)t 



1+ 



(1) 

(2) 
(3) 

(4) 

(5) 



with P* a = (D*°,D t+ ,D* + ) or (B*-,B*°,Bf) being the 
charmed or antibottomed vector mesons andP = (D°,D + ,D + S ) 
or (B~,B°, B®) being the charmed or anti -bottomed pseu- 
doscalar mesons. The trace acts on the gamma matrices. 
The axial-current is defined as A" = \ (<f 3^ - = 

■jd^M + . . . where £ = e f* with M being the exchanged pseu- 
doscalar meson matrix given in Eq. ||6). The vector current 
V is defined as V = \ (^3^ + £5^). In the heavy quark 
limit, the heavy meson velocity is adopted as = (1, 0, 0, 0). 
Fftvip) = d^py - - [p M ,p v ] where p p = ^|p^ with p being 
the exchanged vector meson matrix given in Eq. [7] Expand- 
ing the Lagrangians given in Eqs. dTj3j, we list the specific 
expressions in Eqs. (l8lfT4li. 



At 



V2 V6 
if 



1^0 + jl go 

V2 V6 

K° 



(6) 



P 1 = 



y_ _|_ w 

V2 V2 
P" 

K*- 



K 



.4+ *> K *o 
V5 V2 



(7) 



-Cp>p*m = -i% aflvA xPP*fP*fd v M ba , (8) 

In 

Lp-pm = ~{P b Pl + PliP^d A M b «, (9) 

In 



-Cppv 

-Cp'PV 
■Cp'P'V 



-^2Bg v P b P\vp ba , (10) 
-2 y/UgySew (P b Pf + P'fPl) {d a ff) ba Q D 
yf2£g v P* b -P*Jvp ba 

-i2 -fagvlfpy (drh - d v pX , (12) 



-Cppa = -2g s P b P b <T, 



(13) 
(14) 



In the above, f„ = 132 MeV is the pion decay con- 
stant. The coupling constants g was studied by many theo- 
retical approaches, such as quark model lO and QCD sum 
rule (361 B7F. In our study, we take the experimental re- 
sult of the CLEO Collaboration, g = 0.59 + 0.07 + 0.01, 
which was extracted from the full width of D* + lf38h . For 
the coupling constants relative to the vector meson exchange, 
we adopt the values g v — 5.8 and B = 0.9 which were de- 
termined by the vector meson dominance mechanism, and 
A = 0.56 GeV 1 which was obtained by matching the form 
factor predicted by the effective theory approach with that ob- 
tained by the light cone sum rule and the lattice QCD simu- 
lation 09l l40ll . The coupling constant for the scalar meson 
exchange is g s = g„l(2 Vo~) IToll with g n = 3.73. We take the 
masses of the heavy mesons and the exchanged light mesons 
from PDG Il4lll and summarize them in Table HIT1 



B. The derivation of the effective potentials 

Using the lagrangians given in Eqs. (l8lfT4b. one can easily 
deduce the effective potentials in the momentum space. Tak- 
ing into account the structure effect of the heavy mesons, we 
introduce a monopole form factor 



F(q) 



A 2 - mj 
A 2 -q 2 



(15) 



at each vertex. Here, A is the cutoff parameter and m ex is the 
mass of the exchanged meson. Making fourier transformation 



V(r) 



dq 3 e- itir V(q)F 2 (q) 



(16) 



one can obtain the effective potentials in the coordinate space. 
In Eqs. ( fTSIIBTt , we list the specific expressions of the effec- 
tive subpotentials which are flavor-independent. The effective 



3 



TABLE I: The different channels for the D f *'>D M systems, and similarly, for the B (,) B M . S is the strangeness while / stands for the isospin 
of the systems, "x" means the corresponding state dose not exist due to the symmetry. For simplicity, we adopt the shorthand notations, 
[£»£»*]_ = -4 (DD' - D'D) and [£>£>*]+ = j= (DD* + D'D). 

Channels 



5 


/ 


J 1 ' 


1 


2 


3 


4 


5 6 7 






+ 








X 









1 + 




rnn'i finA 


u u \ o i; 











2+ 


yuu u 2 ) 


n*n*( 3 

1 J i t ^ '-'2/ 










+ 


DD( l S ) 


D'D'( l S ) 


d*d*( 5 a>) 










1 + 
1 

2+ 


r r\ / V' i /3 c \ 
{DD J+( 6 J 

D'D'( 5 S 2 ) 


LOO J+( 

d*d*Cd 2 ) 


D D ( D\) 

D*D*( 5 D 2 ) 










0+ 


DD S CS ) 


D'D' s ( l S ) 


D'D' S ( 5 D ) 






1 


i 

2 


r 

2 + 


DDtfSx) 
D'D' S ( 5 S 2 ) 


DD' S ( 3 D { ) 
D'D' s ( l D 2 ) 


D'D S ( 3 S i) 
£>*£»:( 5 £» 2 ) 




D » D . ( 3 Sl) d-d'^d,) inrtfDi) 






+ 


D s D s ( l S ) 


d:d:( 1 s ) 


d:d:( 5 £» ) 






2 





r 

2 + 


d;d;( 5 s 2 ) 


[D S D' S U 3 D!) 
D' S D' S CD 2 ) 


£>;£>:( 5 £»i) 
£»:£>:( s £» 2 ) 







TABLE II: The different channels for the DW§M systems. "5 1 " is the strangeness of the corresponding system while "/" stands for the isospin 
of the state. 

Channels 

5 / J p 1 2 3 4 5 6 7 8 9 10 

+ £>S('So) D'B'CSo) D'B'( 5 D ) 

1+ DS*( 3 5i) DB'( 3 D t ) D'BCSi) D'B( 3 Di) £)*£*(%) D'B'CD,) D'B'( 5 D { ) 
2 + D'B'( 5 S 2 ) D'B'( [ D 2 ) D'B'( S D 2 ) 

0+ DBOSo) Z)*S*('5o) D'B'( 5 D ) 

1 1 + DB*( 3 Si) DB'( 3 D,) D*B( 3 S{) D'B( 3 D { ) D'B'( 3 S { ) D*B*( 3 Di) D'B'( s D x ) 
2 + d*B*( 5 S 2 ) D'B'( l D 2 ) D'B'( 5 D 2 ) 

0+ Dfl/'So) Dj('So) D'B' S CS ) D' S B'CS ) D'B' S ( 5 D ) D' S B*( 5 D ) 

1 | 1 + D s S*( 3 S0 r>B*( 3 S,) £>:B( 3 5,) ITB,<?Si) £>:B*( 3 5,) D'BlCSi) D^DO DS:( 3 Di) £»:B( 3 £»!) D'B.fDi) 

2 + D;g-( 5 5 2 ) ZTfi^S,) D' s B'( l D 2 ) D'B'^D,) D' s B'eD 2 ) D'Bl( s D 2 ) 

o + D s B„es ) d'M i s ) d:b:( 5 d ) 

2 1+ D^CSi) D s B*( 3 Di) D:B s ( 3 S<) D* s B s ( 3 D0 D' s B' s ( 3 S<) D*B*( 3 Di) D*B*( 5 Di) 

2 + £>:5:( 5 S 2 ) D' s B' s ( l D 2 ) D* S B*( 5 D 2 ) 



potential used in our calculation is the product of the flavor- 
independent subpotentials and the isospin-dependent coeffi- 
cients which are summarized in the AppendixlVIl The flavor- 
independent subpotentials are 

Or) = -«?flb(A,?o,m«r,r), (17) 
V?(r) = ^H (A,q Q ,m v ,r), (18) 



for the process PP — > PP, 



^ Ml -W = 



(19) 



for the scattering process PP — > PP* 



V c (r) = — 

J n 

1 



H 3 (A, q ,m p ,r)T(e\,e\) 



+ -H l (A, q , m p , r)S (e' v e\) 



(20) 
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TABLE III: The masses of the heavy mesons and the exchanged light 
mesons taken from the PDG l4lll . In our study, we keep the isospin 
symmetry. For the isospin multiplet, we use the averaged mass in our 
study. 

Heavy Mesons Mass (MeV) Exchanged Mesons Mass (MeV) 



D ± 


1869.60 


7T ± 


139.57 


D° 


1864.83 


n° 


134.98 


D ,± 


2010.25 


'1 


547.85 


D*° 


2006.96 


P 


775.49 


D: 


1968.47 


CO 


782.65 


D? 


2112.3 


<t> 


1019.46 


B 


5325.1 


<r 


600 


B ± 


5279.17 


A'* 


493.67 


fl° 


5279.50 


K° 


497.61 


«S 


5366.3 


K ,± 


891.66 


B* 


5415.4 


K*° 


895.94 



for the scattering process PP* — » p*p* , and 

2 r 

-j H 3 (A, q , m p , r)T (ie\ x e u ie\ x e 2 ) 
+]-H x (A, q , m p , r)S x e u ie' 4 x e 2 ) 



V h p (r) 



V'lXr) 



= -g 2 H Q (A, q , nto-, r)C {ie\ 



P 2 g 2 y 



X e x ,ie\ X e 2 
H (A, q , m v , r)C (ie| x e\ , ie\ x e 2 ) 
H 3 (A, q , m v , r)T (ic\ x e x ,ie\ x e 2 ) 



(29) 
(30) 



-2A 2 g 2 v 



2 i 
--Hi(A,q ,m v ,r)S [i 



el x e\,ie\ x 62 



(31) 



for the scattering process P*P* — * P*P* . To obtain the effec- 
tive potentials V s {r) for the process P*P — > P*P*, one just 
needs to make the following changes 



V c v (r) = -2A z g z v 



H 3 (A,q Q ,m v ,r)T(e ] v el) 



--Hi (A, tf , m V) r)5 (63, e\) 
for the scattering process PP — > P*.P*, 

V£(r) = -«?#o(A, ?o, m«r,r)S(eJ,ei) 

^2 

-#o(A, 90, m v , r)S (el, e 2 ) 



V d v {r) = 



for the scattering process PP* — > PP*, 



Kir) = 



f 2 

.In 



M 3 (A, qo,m n , r)T(el,e 2 ) 



+ -Mi(A, q , m n , r)S (e\, e 2 ) 



Jn 



H 3 (A, q Q , m„ /r , /K , r)T(e\, e 2 ) 



+ jHi(A,q , m nlnlK , r)S (e\,e 2 ) 



V e Jr) 



o ,2 2 

-2/1 g v 



H-i(A,qo,m v ,r)T(el,e 2 ) 



: Hi(A,qo,m v ,r)S(el,e2) 



3"™' -~ 
for the scattering process PP* — > P*P, 



V f p (r) 



f 2 

J 71 



H 3 (A, q , nip, r)T(e\, ie\ x e 2 ) 



1 



(21) 

(22) 
(23) 



(24) 



(25) 



(26) 



(27) 



+-H X (A, q ,m p ,r)T(e ] v ie\ x e 2 ) 
W(r) = 2^Vv{^3(A,^o,'« l „r) 

x[T(ielxel,e 2 )-T(ielxel,el)] 
1 

+-#i(A, q Q ,m v ,r) 

x[S(iel x e J, e 2 ) - S(ieJ x e\, e J)]} (28) 



6 3 "» f 4' 

62 -> ei, 

6 1 -» e 3- 



164 x e 2 -» ;e 3 X e{, 
ic| X -» icj x e|, 
ie\ Xf 2 -> ie} x ej , 



(32) 



in Eqs.(l27b-d28b. Functions H , H\, H 3 , M\ and M 3 are given 
in the Appendix. Operator C, the generalized tensor operator 
T and spin-spin operator S are defined as 



C(a, b) = ab, 

T(a,b) = a b, 

S(a,b) = a b. 



(33) 
(34) 
(35) 



Due to the large mass gap between the mesons D(D°,D + ), 
£>+, D*(D*°, D* + ) and D* + (similarly, in the bottom sector), it 
is necessary to adopt the nonzero time component of the trans- 
ferred momentum for some scattering processes. We present 
the go s use d in our calculation in the Appendix. Notice that 
mo* - mo > nix leads to the complex potential for the scatter- 
ing process DD* — > D*D, and we take its real part which has 
an oscillation form, see Eq. d24b . 



III. NUMERICAL RESULTS 



Using the potentials given in the subsection III Bl we solve 
the coupled-channel Schrodinger equation and summarize the 
numerical results which include the binding energy (B.E.), 
the system mass (M), the root-mean-square radius (r rms ) and 
the probability of the individual channel (P,) in Tables Hv][V] 

[ylTl rvml Hxb ndlxl 



A. The Numerical Results for Systems with Strangeness 5=0 

For the systems with strangeness S = 0, in order to high- 
light the role of the long-range pion exchange in the formation 
of the loosely bound state, we first give the numerical results 
with the pion-exchange potential alone, which are marked 



5 



with OPE, and then with the heavier eta, sigma, rho and omega 
exchanges as well as the pion exchange, which are marked 
with OBE, see Tables[iVl[V]and|Vl] 



1. D«D<*> 





The state D { * ) D i * ) [I(J p ) = 0(0 + )] is forbidden because the 
present boson system should satisfy the Boson-Einstein statis- 
tic. However, the state D i * ) D { * ) [I(J p ) = 0(1 + )] is very inter- 
esting. Using the long-range pion exchange potential, we ob- 
tain a loosely bound state with a reasonable cutoff. For our 
present D w D w [/( J p ) = 0(1 + )] state, with the cutoff parame- 
ter fixed larger than 1 .05 Ge V, the long-range pion exchange 
is strong enough to form the loosely bound state. If we set 
the cutoff parameter to be 1 .05 Ge V, the binding energy rela- 
tive to the DD* threshold is 1 .24 MeV and the corresponding 
root-mean-square radius is 3 . 1 1 fm which is comparable to the 
size of the deuteron (about 2.0 fm). The dominant channel is 
[DD*]^( 3 Si), with a probability 96.39%. With such a large 
mass gap (about 140 MeV) between the threshold of DD* 
and that of D*D*, the contribution of the state D*D*( 3 S\) is 
2.79%. However, the probability of the D-wave is around 1%. 
When we tune the cutoff to be 1 .20 GeV, the binding energy is 
20.98 MeV and the root-mean-square radius changes into 0.84 
fm. When we use the one-boson-exchange potential, we no- 
tice that the binding becomes deeper. For example, if the cut- 
off is fixed at 1.10 GeV, the binding energy is 4.63 MeV with 
OPE potential. However, it changes into 42.82 GeV with the 
OBE potential for the same cutoff, see Table II VI We also plot 
the potentials in Fig. Q] From the potentials, one can see that 
the heavier rho and omega exchange cancel each other sig- 
nificantly, and the residual interaction is helpful to strengthen 
the binding. The present meson exchange approach strongly 
suggests D ( * ) D ( * ) [I(J P ) = 0(1 + )] to be a molecule. 

We should mention that in the calculation of the X(3872) 
one also obtained a bound DD* state with quantum numbers 
I(J PC ) = 0(1 ++ ) using the OPE potential EH- One may 
be confused since the difference between the potential of the 
DD* system and that of the DD* system is the G-parity of the 
exchanged meson while the pion has an odd G-parity. Actu- 
ally, the iso-singlet DD* system has two C-parity states, one 
with even C-parity (C = +) and the other with odd C-parity 
(C = -). And, the interaction of our present DD* system re- 
lates to that of the odd C-parity but not the even C-parity DD* 
state via the G-parity rule. 

We obtain no binding solutions for the state 
DW£)W[/(y p ) = 0(2+)] even if we tune the cutoff pa- 
rameter as high as 3.0 GeV. Therefore, the existence of the 
state D^ D ( -*\l{J p ) = 0(2 + )] as a molecule is ruled out within 
the present meson exchange approach. 

For the state D^* ) D i * ) [I(J p ) = 1(0 + )], with the cutoff less 
than 3.0 GeV, the long-range pion exchange is not sufficient to 
form the bound state. However, when we add the heavier eta, 
sigma, rho and omega exchanges and tune the cutoff to be 2.64 
MeV, a bound state with mass 3730.17 MeV appears. The 
binding energy relative to the DD threshold is 2.64 MeV and 
the corresponding root-mean-square radius is 1.38 fm. The 



(Vn) 



(V22) 




(V33) 




(V 12 ) 





(Vm) 



(V44) 




(V13) 




(V23) 




(V 34 ) 



FIG. 1: Color online. The effective potentials for the state 
D <t) D < * ) [I(J p ) = 0(1 + )] with cutoff = 1.00 GeV. 



channel DD( Sq) with a probability of 91.47% dominates this 
state. The probability of the channel D*D*( 5 Dq) is very small, 
only 0.07%. 

Just as the J p = + case, with the pion-exchange po- 
tential alone we obtain no binding solutions for the state 
D ( * ) D ( * ) [I(J P ) = 1(1 + )] with the cutoff parameter less than 
3.0 GeV. When we use the OBE potential and tune the cutoff 
to be 2.48 GeV, we obtain a bound D ( * ) D ( * ) [I(J P ) = 1(1 + )] 
state with mass 3875.58 MeV. The binding energy relative to 
the DD* threshold is 0.27 MeV and the corresponding root- 
mean-square radius is 5.81 fm. The dominant channel is 
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[DZ)*] + ( 3 5i), with a probability of 99.92%. 

For the state D w D w [/(./ p ) = 1(2 + )], it is necessary to 
mention that there are five channels DD( l D2), [DD*] + ( 3 D2), 
D*D*( 5 S 2 ), D*D*( 1 D 2 ) and D*D*( 5 D 2 ), with the quantum 
numbers I(J P ) = 1(2 + ). If we consider all the five channels, 
with the OBE potential we obtain a bound state with the cutoff 
parameter fixed to be 2.84 GeV. The binding energy relative to 
the DD threshold is 12.93 MeV. Surprisingly, the correspond- 
ing root-mean-square radius is as small as 0.22 fm. The domi- 
nant channel is D*D*( 5 5 , 2), with a probability of 99.5%. How- 
ever, the probability of the channel DD( l D2) is so small as 
0.04%, which tells us that this is not a loosely bound DD state 
but a deeply bound D*D* state. With so tight a bound state, 
the present meson-exchange model dose not work. Therefore, 
we omit the channels DD( i D 2 ) and [DD*] + ( 3 D2) and keep the 
three D*D* channels. With the pion-exchange potential we 
fail to obtain a bound state with the cutoff parameter less than 
3.0 GeV. However, when we use the OBE potential and tune 
the cutoff to be 2.48 GeV, we obtain a bound state with mass 
4014.29 MeV. The binding energy is 2.95 MeV and the cor- 
responding root-mean-square radius is 1.61 fm. The channel 
D*D*( 5 Sz) with a probability of 99.93% dominates this state. 

For the states D W D ( * ) [/(7 P ) = 1(0 + ), 1(1 + ), 1(2 + )], with the 
present meson exchange approach we can not draw a definite 
conclusion about whether they are molecules or not. Further 
detailed study with other approaches are needed. 



2. 

With the heavy quark flavor symmetry, the potentials for the 
gMjjO) system are similar to those for the DD system. The 
main difference between the two systems is that the reduced 
mass of the B^S" system is much larger than that of the 
£)(*)£)(*) S y S tem. We summarize our numerical results of the 
SMgOO system in Table [V] 

Similar to the charmed state D w D w [/(/ p ) = 0(1 + )], the 
bottomed state B { * ) B {t) [I(J p ) = 0(1 + )] is also very interest- 
ing. The long-range pion exchange is strong enough to form 
the loosely bound B i * ) B { * ) [I{J p ) = 0(1 + )] state with the cutoff 
larger than 1 .25 GeV. If we tune the cutoff to be 1.30 GeV, the 
the binding energy is 1.14 MeV and the corresponding root- 
mean-square radius is 2.25 fm. The channel [BB*]-( 3 S i) with 
a probability of 91.83% dominates this state. The probability 
of the D-wave is 6.96%, see Table [Vl When we use the OBE 
potential, the binding becomes tighter as expected, which is 
similar to its charmed partner D {t) D {t) [I(J p ) = 0(1 + )]. With 
the present meson exchange approach we also suggest that the 
state B ( * ) B i * ) [I(J p ) = 0(1 + )] is a good candidate of molecule. 

Different from its charmed partner, B^*>B ( -* ) [I(J P ) = 0(2 + )] 
can form a bound state with the pion-exchange potential if the 
cutoff is tuned larger than 2.88 GeV. When we set the cutoff 
to be 2.88 GeV, the binding energy is 2.75 MeV, and corre- 
spondingly, the root-mean-square radius is 0.72 fm. When we 
use the OBE potential, we obtain the binding solutions with a 
smaller but more reasonable cutoff. Unfortunately, the bind- 
ing solutions depend very sensitively on the cutoff parameter. 
When tune the cutoff from 1 .66 GeV to 1 .72 GeV, the binding 



energy changes from 8.30 MeV to 73.89 MeV. The numerical 
results with the present meson exchange approach suggests 
that the states B (t) B ( *\l(J p ) = 0(2 + )] can also be a candidate 
of molecule, but not a good one because of the strong depen- 
dence of the binding solutions on the cutoff. 

The pion-exchange alone is also sufficient to form the 
loosely bound B ( * ) B { * ) [I(J P ) = 1(0 + )] state with the cutoff 
larger than 1 .70 GeV. When we tune the cutoff from 1 .70 GeV 
to 1 .90 GeV, the binding energy increases from 1 .05 MeV to 
1 1 .29 MeV, and correspondingly, the root-mean-square radius 
decreases from 2.07 fm to 0.75 fm. The dominant channel 
is BB( l S ), with a probability of 95.35% ~ 90.29%. And, 
the probability of the channel D*D*( 5 D ) is 1.79% ~ 4.69%. 
Actually, two pseudoscalar D-mesons can not interact with 
each other via exchanging a pion. Therefore, the binding so- 
lutions totally come from the coupled-channel effect, just as 
in the AgAg case l45l l46ll . When we add the contributions of 
the heavier eta, sigma, rho and omega exchanges, the results 
change little, which implies that the eta, sigma, rho and omega 
exchanges cancel with each other significantly. The existence 
of the state B ( * ) B i * ) [I(J p ) = 1(0 + )] as a molecule is favored by 
the present numerical results. 

The state B i * ) B ( * ) [I(J p ) = 1(1 + )] is also an interesting one. 
When we tune the cutoff between 1 .40 GeV and 1 .70 GeV, 
we also obtain a loosely bound state with the OPE potential. 
The binding energy is 0.83 ~ 1 1.80 MeV and the correspond- 
ing root-mean-square radius is 2.36 ~ 0.80 fm. The domi- 
nant channel is [BB*]+( 3 Si), with a probability of 96.59% ~ 
90.63%. Different from the isospin singlet case, when we use 
the OBE potential the binding becomes shallower. For exam- 
ple, with the OPE potential the binding energy is 1 1 .80 MeV 
if the cutoff is set to be 1 .70 MeV while with the OBE po- 
tential it is 8.29 MeV for the same cutoff, see Table M The 
present numerical results suggest that B ( * ) B ( * ) [I(J P ) = 1(1 + )] 
can also be a good candidate of molecule. 

With the same reason as for the charmed case, we only keep 
the three channels B*B*( 5 S 2 ), B*B*( 1 D 2 ) and B*B*( 5 D 2 ) for 
the state B { * ) B ( * ) [I(J P ) = 0(2 + )]. with the pion-exchange po- 
tential, when we tune the cutoff between 1.80 GeV and 2.10 
GeV, we obtain a loosely bound state with binding energy 
2.25 ~ 13.35 MeV and root-mean-square radius 1.48 ~ 0.71 
fm. The dominant channel is B*B*( 5 5 , 2), with a probability 
of 95.73% ~ 93.66%. When we use the OBE potential and 
tune the cutoff from 1 .70 GeV to 1 .76 GeV, the binding en- 
ergy changes into 7.88 ~ 28.01 MeV, and correspondingly, 
the root-mean-square radius changes into 0.59 ~ 0.36 fm. The 
existence of B (,) B^[I(J P ) = 1(2 + )] as a molecule is also fa- 
vored by our present meson exchange model. 



3. 

The small binding energy and large root-mean-square 
radius with reasonable cutoff parameter makes the state 
D { * ) B { *\I(J P ) = 0(0 + )] very interesting. With the OPE po- 
tential, when we fix the cutoff between 1 .08 GeV and 1 .20 
GeV, we obtain a loosely bound state with binding energy 
0.42 ~ 31.35 MeV and root-mean-square radius 4.23 ~ 0.55 
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TABLE IV: The numerical results for the D { ' ] D { ' ] system with strangeness 


5=0. "x" means the corresponds 


ig state dose not exit due to 


symmetry while "-" means there dose not exist binding energy with the cutoff parameter less than 3.0 GeV. The binding energies for the states 


D M D M [I(JP) 


= 0(l + )]and D { 


,} D M [I(J P ) = 1(1+)] are relative to the threshold of DD' while that of the state D i ' ) D { ' ) [I(J p ) = 1(0 + )] is relative 


to the DD threshold. 


















£)(*)£)(*) 








/ J p 




OPE 






OBE 




o+ 




X 






X 






A(GeV) 


1 IK 1 If! 1 IS 


1 9n n os 

l.^u \J.yJ 


1 on 


1 .UJ 


1.1U 




B.E. (MeV) 


1.24 4.63 11.02 


20.98 0.47 


5.44 


18.72 


42.82 




M (MeV) 


3874.61 3871.22 3864.83 


3854.87 3875.38 


3870.41 


3857.13 


3833.03 


r 


r rml (fm) 


3.11 1.68 1.12 


0.84 4.46 


1.58 


0.91 


0.64 




Pj(%) 


96.39 92.71 88.22 


83.34 97.97 


92.94 


85.64 


77.88 




Pz(%) 


0.73 0.72 0.57 


0.42 0.58 


0.55 


0.32 


0.15 




P 3 (%) 


2.79 6.45 11.07 


16.11 1.41 


6.42 


13.97 


21.91 




Pa{%) 


0.08 0.13 0.14 


0.13 0.04 


0.09 


0.08 


0.05 


2+ ........ 




A(GeV) 




2.64 


2.66 


2.68 


2.70 




B.E. (MeV) 




4.29 


12.63 


23.30 


35.75 


+ 


M(MeV) 




3730.17 


3721.83 


3711.16 


3698.71 


r rms (fm) 




1.38 


0.79 


0.58 


0.48 




Pi(%) 




91.47 


88.26 


86.25 


84.80 




Pi(%) 




8.46 


11.64 


13.63 


15.07 




P 3 (%) 




0.07 


0.10 


0.12 


0.13 




A(GeV) 




2.48 


2.50 


2.52 


2.54 




B.E. (MeV) 




0.27 


3.82 


9.79 


17.39 


1 1+ 


M(MeV) 




3875.58 


3872.03 


3866.06 


3858.46 


r rm s(fm) 




5.81 


1.47 


0.91 


0.69 




Pz(.%) 




99.92 


99.90 


99.90 


99.90 




P 3 (.%) 




0.06 


0.06 


0.04 


0.04 




Pa(%) 




0.02 


0.05 


0.06 


0.06 




A(GeV) 




2.48 


2.50 


2.52 


2.54 




B.E. (MeV) 




2.95 


8.86 


16.51 


25.54 


2 + 


M(MeV) 




4014.29 


4008.38 


4000.73 


3991.70 




r rms (fm) 




1.61 


0.92 


0.68 


0.56 




Pi(.%) 




99.93 


99.94 


99.95 


99.95 




P%{%) 




0.01 


0.01 


0.01 


0.00 




P 3 (%) 




0.06 


0.05 


0.04 


0.04 



fm. The dominant channel is DB{ 1 Sq), with a probability of 
95.32% ~ 60.03%. The probability of the channel D*B*( 5 D ) 
is very small as expected, see TablelVIl When we add the con- 
tributions of the heavier eta, sigma, rho and omega exchanges, 
the binding energy changes by tens of MeV. With the present 
numerical results of the state D ( * ) B ( * ) [I(J P ) = 0(0 + )], we pre- 
dict it to be a molecular state. 

When the cutoff is tuned larger than 2.05 GeV, the long- 
range pion exchange is strong enough to form the loosely 
bound D { * ) B { * ) [I(J P ) = 0(1 + )] state. If we set the cutoff be- 
tween 2.05 GeV and 2.20 GeV, the binding energy relative to 
the DB* threshold is 1.21 ~ 6.30 MeV while the root-mean- 
square radius is 2.75 ~ 1.37 fm. The dominant channel is 
DB*( 3 Si), with a probability of 96.81% ~ 93.31%. When we 
add the contributions of the heavier eta, sigma, rho and omega 
exchanges, we obtain a loosely bound D ( * ) B ( * ) [/(7 F ) = 0(1 + )] 



state with a reasonable cutoff 1.65 ~ 1.80 GeV. If we set 
the cutoff parameter to be 1.70 GeV, the binding energy is 
2.83 MeV and the root-mean-square radius is 1.89 fm which 
is comparable to the size of the deuteron (about 2.0 fm). 
The channel DB*( 3 Si) with a probability of 95.20% dom- 
inates this state. However, the contribution of the D-wave 
is small, less than 5.0%. It seems that the possibility of 
£)(*)/}(*) [I(J P ) = 0(1 + )] as a molecule is favored by the present 
meson exchange approach. 

With the same reason as for the state B ( ^B (t, [I(J p ) - 
0(2+)], we omit the channels DB{ l D 2 ), DB*( 3 D 2 ) and 
D*B( 3 D 2 ) for the state D ( * ) B<-* ) [I(J P ) = 0(2 + )]. Since the 
amplitudes of the channel D*B*( 1 D2) scattering into the other 
channels are zero, we also omit this channel in our study. With 
the OPE potential, when we set the cutoff between 2.10 GeV 
and 2.30 GeV, we obtain a bound D^B^[I(J P ) = 0(2 + )] 
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TABLE V: The numerical results the flW/fW systems with strangeness 5=0. "x" means the corresponding state dose not exist due to the 
symmetry. The binding energies of the states B M B M [l(J p ) = 0(1 + )], B ( * ) B W [7(7 P ) = 0(2+)] and B {,) B {, \I{J P ) = 1(1 + )] are relative to the BB* 
threshold while that of the state B f ' ) B i ' ) [I(J p ) = 1(0 + )] is relative to the BB threshold. 

/ J p OPE OBE 

+ x x 





A(GeV) 


1.25 


1.30 


1.40 


1.50 


1.10 


1.15 


1.20 


1.25 




B.E. (MeV) 


0.36 


1.14 


4.10 


9.44 


0.19 


1.56 


5.20 


13.71 




M (MeV) 


10604.08 


10603.30 


10600.34 


10595.00 


10604.25 


10602.88 


10599.24 


10590.73 


r 


r rms (fm) 


3.60 


2.25 


1.40 


1.04 


4.86 


1.99 


1.26 


0.86 




Pi(%) 


94.29 


91.83 


89.03 


83.18 


96.66 


94.17 


92.34 


86.54 




Pi(%) 


4.16 


5.66 


7.01 


7.10 


2.61 


4.04 


3.54 


1.96 




PA%) 


0.77 


1.20 


1.74 


1.89 


0.36 


0.97 


3.14 


10.79 




P 4 (%) 


0.78 


1.30 


2.21 


2.83 


0.37 


0.82 


0.98 


0.72 




A(GeV) 


2.88 


2.90 


2.94 


2.96 


1.66 


1.68 


1.70 


1.72 




B.E. (MeV) 


2.75 


5.18 


10.38 


13.15 


8.30 


28.01 


49.89 


73.89 


2* 


r rms fm 


0.72 


0.68 


0.63 


0.61 


0.39 


0.35 


0.32 


0.31 




M(MeV) 


10601.69 


10599.26 


10594.06 


10591.29 


10596.14 


10576.43 


10554.55 


10530.55 




Pi(%) 


61.15 


60.47 


59.43 


59.01 


57.82 


55.86 


54.69 


53.88 




Pi(%) 


38.85 


39.53 


40.57 


40.99 


42.18 


44.14 


45.31 


46.12 




A(GeV) 


1.70 


1.75 


1.80 


1.90 


1.74 


1.76 


1.78 


1.80 




B.E. (MeV) 


1.05 


2.53 


4.70 


11.29 


2.24 


8.19 


15.67 


24.33 




M(MeV) 


10557.63 


10556.15 


10553.98 


10547.39 


10556.44 


10550.49 


10543.01 


10534.35 


+ 


r rml (fm) 


2.07 


1.40 


1.07 


0.75 


1.03 


0.55 


0.42 


0.36 




Pi{%) 


95.35 


92.82 


90.29 


85.46 


87.27 


83.59 


81.81 


80.66 




Pi(%) 


2.02 


3.20 


2.42 


6.84 


12.54 


16.26 


18.07 


19.25 




P 3 (%) 


2.63 


3.98 


5.28 


7.70 


0.19 


0.15 


0.11 


0.08 




A(GeV) 


1.40 


1.50 


1.60 


1.70 


1.66 


1.68 


1.70 


1.72 




B.E. (MeV) 


0.83 


1.89 


6.46 


11.80 


0.55 


3.72 


8.29 


13.92 




M(MeV) 


10603.61 


10601.55 


10597.98 


10592.64 


10603.89 


10600.72 


10596.15 


10590.52 


r 


r rms (fm) 


2.36 


1.38 


1.00 


0.80 


2.26 


0.85 


0.59 


0.48 




Pi(%) 


96.59 


94.40 


92.43 


90.63 


99.24 


99.34 


99.51 


99.64 




Pi(%) 


1.62 


2.35 


2.88 


3.29 


0.41 


0.30 


0.20 


0.14 




P 3 (%) 


1.79 


3.25 


4.69 


6.08 


0.34 


0.36 


0.29 


0.22 




A(GeV) 


1.80 


1.90 


2.00 


2.10 


1.70 


1.72 


1.74 


1.76 




B.E. (MeV) 


2.25 


4.74 


8.38 


13.35 


7.88 


13.63 


20.36 


28.01 




M(MeV) 


10647.95 


10645.46 


10641.82 


10636.85 


10642.32 


10636.57 


10629.84 


10622.19 




' rms 


1.48 


1.08 


0.85 


0.71 


0.59 


0.47 


0.41 


0.36 




Pi(%) 


95.73 


94.86 


94.19 


93.66 


99.70 


99.80 


99.86 


99.91 




p 2 (%) 


0.72 


0.87 


0.97 


1.06 


0.05 


0.03 


0.02 


0.01 




p 3 (.%) 


3.55 


4.27 


4.83 


5.28 


0.25 


0.17 


0.12 


0.08 



state, with binding energy 0.80 ~ 3.83 fm and root-mean- 
square radius 3.35 ~ 1.73 fm. The channel D*B*( 5 S2) pro- 
vides a dominant contribution, 94.14% ~ 89.27%. When we 
use the OBE potential, the binding energy is 0.63 ~ 1 .79 MeV 
and the root-mean-square radius is 3.70 ~ 2.43 fm for the cut- 
off between 1.90 GeV and 2.10 GeV. Such a loosely bound 
state with weak dependence of the binding solutions on the 
cutoff parameter is particularly interesting. And, it seems that 
our present meson exchange model favors the existence of 
D { * } B[I(J P ) = 0(2+)] as a molecular state. 

For the state D < * ) B < * ) [I(J P ) = 1(0 + )], if we tune the cutoff 



parameter larger than 2.60 GeV, the OPE potential is sufficient 
to form the D ( * ) B ( * ) [I(J P ) = 1(0 + )] bound state. When we 
fix the cutoff parameter between 2.60 GeV and 2.90 Gev, the 
binding energy is 0.14 ~ 17.53 fm and the corresponding root- 
mean-square radius is 7.05 ~ 0.75 fm. The dominant channel 
is DB( l S ), with a probability of 98.68% ~ 87.27%. When 
we use the OBE potential, we obtain binding solutions with a 
cutoff parameter larger than 2.22 GeV, see TablelVIl However, 
with such a large cutolf parameter, we can not draw a definite 
conclusion about whether the state D W B W [7(/ P ) = 1(0 + )] is 
a molecule or not. And, further detailed study with other ap- 
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proach will be helpful to settle this issue. 

For the state D w 5 w [/(./ p ) = 1(1 + )], when we tune the cut- 
off parameter larger than 2.55 GeV, we obtain binding solu- 
tions with the OPE potential. If we set the cutoff parameter to 
be 2.60 GeV, the binding energy relative to the DB* threshold 
is 2.43 MeV and the corresponding root-mean-square radius 
is 1.83 fm. However, when we add the heavier eta, sigma, rho 
and omega exchanges, we obtain no binding solutions with 
the cutoff parameter less than 3.0 GeV. It seems that the pos- 
sibility of the state D { * ) B i * ) [I(J p ) = 1(1 + )] as a molecule is 
ruled out by the present meson exchange model. 

For the state = 1(2 + )], when we tune the cut- 

off parameter as large as 2.80 GeV, we obtain binding solu- 
tions with the OPE potential. If we set the cutoff parameter 
to be 2.90 GeV, the binding energy is 2.00 MeV and the cor- 
responding root-mean-square radius is 1.95 fm. The domi- 
nant channel is D*B*( 5 S 2 ), with a probability of 97.56%. The 
probability of the D-wave is 2.44%. When we use the OBE 
potential, we obtain binding solutions with a smaller cutoff, 
see TablelVIl If we tune the cutoff to be 2.10 GeV, the binding 
energy is 0.44 MeV. Similar to the I(J P ) = 1(0 + ) case, we can 
not draw a definite conclusion with such a large cutoff. And, 
further detailed study with other approach is needed. 



B. The Results for The Systems with Strangeness S = 1 

For the systems with strangeness S = 1 , there dose not exist 
the long-range pion exchange, but there are heavier eta, sigma, 
K and K* exchanges. We summarize the our numerical results 
in Table|Vn]for (DD) S and (BB) S and in Table[VlII]for (OS),. 



I (D { ' } D {,> ) S and (S W S W ) S 

When we fix the cutoff parameter between 2.70 GeV and 
2.76 GeV, we obtain a bound state of (D ( ' ) D ( ' ) ) S [J P = + ] 
with mass between 3832.06 MeV and 3802.15 MeV and a 
large root-mean-square radius 1.53 ~ 0.50 fm. The dom- 
inant channel is DD s ( l So) with a probability of 92.85% ~ 
85.82% and the second largest channel D*D*( l So) contributes 
7.10% ~ 14.07%. The contribution of the channel D*D* S ( 5 D ) 
is very small, only 0.05% ~ 0.11%. Unfortunately, for such 
a loosely bound state, we can not draw a definite conclusion 
with the present meson exchange model about whether it ex- 
ists as a molecule or not. Now the cutoff parameter is as large 
as 2.70 GeV which is twice as large as that for the deuteron 
(1.2 ~ 1.5 GeV l43lo while the binding solutions depend sen- 
sitively on the cutoff parameter. 

In the corresponding bottomed sector, we also obtain a 
bound state of (B ( * ) B ( * ) ) J [7 F = + ] with a smaller and more 
reasonable cutoff, 1.82 ~ 1.88 GeV. This can be easily under- 
stood given that the mass of the bottomed mesons are much 
larger than those of the charmed mesons and the effective po- 
tentials for the two systems are similar. When we tune the 
cutoff to be 1.82 GeV, the binding energy is 0.56 GeV and the 
root-mean-square is 2.28 fm which is comparable to that of 
the deuteron (about 2.0 fm), see Table [VTTl With the present 



meson approach, this state can be a candidate of molecular 
state although not a good one. 

The shallow binding solutions and reasonable cutoff param- 
eter (1.44 ~ 1.50 GeV for (Z) W D W ). S [7 F = 1 + ] and 1.10 ~ 
1.16 GeV for (B M B ( -* ) ) S [J P = 1 + ]), make these two states par- 
ticularly interesting. For the charmed mesons, when we set 
the cutoff to be 1.44 GeV, the binding energy is 5.43 MeV 
and the root-mean-square radius is 1 .36 fm, see Table I VII I 
The closeness of the thresholds for DD* and D*D S makes the 
channels DD* s ( 3 Si) and D*D s ( 3 Si) provide the comparable 
and main contributions, 45.80% for DD*( 3 S } ) and 51.31% 
for D*D s ( 3 Si). The D-wave channel D*D*( 3 Di) and D*D* 
provide almost vanishing contributions, 0.02% for D*D*{3D\) 
and 0.00% for D*D* s { 5 D l ). This is because of the large mass 
gap between D*D* and D*D S and the strong repulsive inter- 
action coming from the centrifugal potential of the D-wave. 
In the bottomed sector, when we fix the cutoff parameter to 
be 1.10 GeV, the binding energy is 0.67 MeV and the root- 
mean-square radius is 2.19 fm. Just as in the charmed case, 
the largest contribution comes from the channels BB*( 3 S\) 
(24.52%) and B*B s ( 3 5i) (72.97%). The probabilities of the 
two channels B*B*( 3 D\) and B*B*( 5 Di) are negligible, either, 
with the same reason as for the charmed case. With the present 
meson exchange approach, it is natural to interpret these two 
states as molecular states. And, we expect these two states 
will be studied in the further experiment. 

For the state (D ( * } D ( * } ) S [J P = 2 + ], we obtain a bound state 
with binding energy 1.54 ~ 22.29 MeV when we tune the 
cutoff to be 2.54 ~ 2.60 GeV. The root-mean square radius is 
2.28 ~ 0.60 fm. The dominant channel is D*D*( 5 S2), with a 
probability of 99.98%. However, the contributions of the D- 
wave are negligible. As mentioned before, for such a bound 
state with a cutoff parameter as large as 2.54 ~ 2.60 GeV, 
we can not draw a definite conclusion with the present meson 
exchange model. 

For the state (B { * ) B ( * ) ) S [J P = 2 + ], if we fix the cutoff pa- 
rameter between 1.76 GeV and 1.82 GeV, the binding en- 
ergy is 0.92 ~ 16.34 MeV and the root-mean-square radius 
is 1.70 ~ 0.44 fm. The dominant channel is B*B*( 5 S2) with 
a probability of 99.76% ~ 99.85%. The present model also 
favors (but not strongly) that this state exist as a molecule. 



2. (D w 5 ( *») s 

For the state (D i * ) B ( * ) ) s [J P = + ], if we set the cutoff 
parameter to be 1.28 GeV, the binding energy is 6.72 MeV 
and correspondingly, the root-mean-square radius is 0.92 fm. 
The probability of the channel DB s CSo) is 50.10% while that 
of the channel D s B( l So) is 25.66%. However, if the bind- 
ing energy is 68.73 MeV, these two channels provide com- 
parable contributions, 25.30% for DB^So) and 23.02% for 
D s B( 1 Sq). Given that the mass gap between DB S and D S B is 
around 14.5 MeV, for the binding energy comparable to this 
value, the mass gap plays an important role in the formation of 
the loosely bound state. However, if the binding energy is as 
large as 68.73 MeV, which is much larger than the mass gap, 
the important effect of the mass gap is gone, which is similar 
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TABLE VI: The numerical results for the £) W B W system with strangeness 5=0. "-" means we obtain no binding solutions for the correspond- 
ing state with the cutoff parameter less than 3.0 GeV. The binding energies of the states D i ' ) B M [I(J p ) = 0(0 + )] and D M B i * ) [l(J p ) = 1(0 + )] are 
relative to the DB' threshold while those of the states D M B M [I(J P ) - 0(1 + )] and D M B ( * ) [I(J P ) = 1(1 + )] are relative to the DB threshold. 

/ J p OPE OBE 









A(GeV) 


1.08 


1.12 


1.16 


1.20 


1.00 


1.02 


1.04 


1.06 




B.E.(MeV) 


0.42 


5.13 


15.54 


31.35 


1.92 


6.84 


15.25 


27.16 


() + 


M(MeV) 


7146.15 


7141.15 


7131.03 


7115.22 


7144.65 


7139.73 


7131.32 


7119.41 


r rms (fm) 


4.23 


1.27 


0.75 


0.55 


2.07 


1.14 


0.78 


0.61 




Pi(%) 


95.32 


82.45 


69.99 


60.03 


91.21 


82.14 


72.93 


64.72 




Pz(%) 


4.39 


17.00 


29.51 


39.59 


8.47 


17.49 


26.75 


35.03 




Pi(%) 


0.28 


0.55 


0.50 


0.38 


0.32 


0.37 


0.32 


0.25 




A(GeV) 


2.05 


2.10 


2.15 


2.20 


1.65 


1.70 


1.75 


1.80 




B.E. (MeV) 


1.21 


2.44 


4.13 


6.30 


0.44 


2.83 


10.75 


30.86 




M (MeV) 


7191.12 


7189.89 


7188.20 


7186.03 


7191.89 


7189.50 


7181.58 


7161.47 




r rms (fm) 


2.75 


2.02 


1.63 


1.37 


4.37 


1.89 


1.05 


0.65 


1+ 


Pi(%) 


96.81 


95.59 


94.42 


93.31 


98.72 


95.20 


84.46 


65.15 




Pi(%) 


0.00 


0.01 


0.01 


0.01 


0.00 


0.00 


0.01 


0.01 




PA%) 


0.22 


0.30 


0.38 


0.46 


0.26 


1.49 


5.96 


13.96 




Pa{%) 


1.10 


1.51 


1.90 


2.26 


0.28 


0.53 


0.58 


0.34 




P S (%) 


0.29 


0.40 


0.50 


0.59 


0.39 


2.08 


8.23 


20.08 




P 6 (%) 


0.38 


0.52 


0.67 


0.80 


0.09 


0.20 


0.24 


0.16 




P 7 (%) 


1.20 


1.67 


2.13 


2.56 


0.26 


0.50 


0.52 


0.29 




A(GeV) 


2.10 


2.15 


2.20 


2.30 


1.90 


1.95 


2.00 


2.10 




B.E. (MeV) 


0.80 


1.31 


1.97 


3.83 


0.63 


0.87 


1.14 


1.79 




M(MeV) 


7332.92 


7332.41 


7331.75 


7329.89 


7333.09 


7332.85 


7332.58 


7331.93 


2 + 


r rms (fm) 


3.35 


2.71 


2.27 


1.73 


3.79 


3.29 


2.93 


2.43 




Pi{%) 


94.14 


92.83 


91.58 


89.27 


96.83 


96.41 


96.03 


95.36 




Pz(%) 


0.25 


0.29 


0.33 


0.38 


0.20 


0.22 


0.24 


0.28 




P 3 (%) 


5.60 


6.88 


8.09 


10.35 


2.97 


3.36 


3.73 


4.36 




A(GeV) 


2.60 


2.70 


2.80 


2.90 


2.22 


2.24 


2.26 


2.28 




B.E. (MeV) 


0.14 


2.56 


8.27 


17.53 


0.62 


6.64 


15.45 


26.01 


() + 


M(MeV) 


7146.43 


7144 01 


7138.30 


7129.04 


7145.95 


7139.93 


7131.12 


7120.56 


r rms (im) 


7.05 


1.78 


1.04 


0.75 


3.05 


0.87 


0.57 


0.45 




Pi(%) 


98.68 


94.88 


90.98 


87.27 


94.14 


87.90 


85.29 


83.65 




P 2 (%) 


0.54 


2.39 


4.36 


6.32 


5.83 


12.05 


14.65 


16.29 




p 3 (%) 


0.66 


2.73 


4.66 


6.41 


0.03 


0.05 


0.06 


0.06 




A(GeV) 


2.55 


2.60 


2.65 


2.70 












B.E. (MeV) 


0.94 


2.43 


4.67 


7.67 


_ 


_ 


_ 


_ 




M (MeV) 


7191.39 


7189.90 


7187.66 


7184.66 


_ 


_ 


_ 


_ 




r rms (fm) 


2.87 


1.83 


1.35 


1.08 


_ 


_ 


_ 


_ 


1+ 


Pi(%) 
P 2 (%) 
P 3 (%) 
P 4 (%) 
Ps(%) 
Pe(%) 
P 7 (%) 


96.13 
0.00 
0.99 
0.70 
0.92 
0.26 
0.99 


93.82 

0.01 

1.64 

1.09 

1.48 

0.41 

1.56 


91.53 
0.01 
2.33 
1.45 
2.03 
0.54 
2.12 


89.31 
0.01 
3.03 
1.77 
2.57 
0.66 
2.66 


_ 
_ 
_ 
_ 


_ 
_ 
_ 
_ 


_ 
_ 
_ 
_ 


_ 
_ 
_ 
_ 




A(GeV) 


2.80 


2.85 


2.90 


2.95 


2.10 


2.12 


2.14 


2.16 




B.E. (MeV) 


0.49 


1.11 


2.00 


3.17 


0.44 


4.64 


10.86 


18.45 




M(MeV) 


7333.23 


7332.61 


7331.72 


7330.55 


7333.28 


7329.08 


7322.86 


7315.27 


2* 


r rms (fm) 


3.80 


2.57 


1.95 


1.58 


3.51 


1.03 


0.68 


0.54 




P a (%) 


98.57 


98.03 


97.56 


97.15 


99.86 


99.89 


99.93 


99.95 




P 2 {%) 


0.24 


0.33 


0.40 


0.47 


0.02 


0.02 


0.01 


0.01 




P 3 (%) 


1.20 


1.65 


2.04 


2.38 


0.12 


0.09 


0.06 


0.04 
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TABLE VII: The numerical results for the D w / B w B M systems with strangeness S = 1. The binding energy of the state 
(D M D ( *>)J(B M B M ) S [J P _ +] is re i at i ve t0 trie DD S threshold while that of the state (D i ' ) D w ) s l(B M &' ) ) s [J P = 1 + ] is relative to the 
DWJB'B, threshold. 





A(GeV) 


2.70 


2.72 


2.74 


2.76 


1.82 


1.84 


1.86 


1.88 




B.E.(MeV) 


3.66 


11.31 


21.47 


33.57 


0.56 


5.27 


12.08 


20.32 


+ 


M(MeV) 


3832.06 


3824.41 


3814.25 


3802.15 


10645.08 


10640.37 


10633.56 


10625.32 


r rms (im) 


1.53 


0.85 


0.62 


0.50 


2.28 


0.71 


0.48 


0.39 




Pi(%) 


92.85 


89.58 


87.42 


85.82 


92.78 


86.34 


83.69 


82.06 




P 2 (%) 


7.10 


10.34 


12.49 


14.07 


7.12 


13.54 


16.21 


17.86 




P3(%) 


0.05 


0.08 


0.09 


0.11 


0.10 


0.12 


0.10 


0.07 




A(GeV) 


1.44 


1.46 


1.48 


1.50 


1.10 


1.12 


1.14 


1.16 




B.E. (MeV) 


5.43 


10.19 


16.31 


23.84 


0.67 


3.21 


7.17 


12.52 




M (MeV) 


3971.68 


3966.92 


3960.80 


3953.27 


10690.73 


10688.19 


10684.23 


10678.88 




r rms (fin) 


1.36 


1.05 


0.87 


0.74 


2.19 


1.10 


0.81 


0.66 


1 + 


Pj(%) 


45.80 


47.33 


47.81 


47.85 


24.53 


36.91 


41.11 


42.75 




P 2 (%) 


0.14 


0.15 


0.15 


0.14 


0.22 


0.32 


0.34 


0.33 




P 3 (%) 


51.31 


48.83 


47.40 


46.41 


72.97 


58.66 


52.92 


49.87 




P 4 (%) 


0.12 


0.13 


0.13 


0.13 


0.36 


0.42 


0.40 


0.36 




p 5 (%) 


2.61 


3.54 


4.49 


5.44 


1.83 


3.55 


5.09 


6.54 




p 6 (%) 


0.02 


0.02 


0.02 


0.02 


0.09 


0.13 


0.14 


0.15 




p 7 (%) 


0.00 


0.00 


0.00 


0.00 


0.01 


0.00 


0.00 


0.00 




A(GeV) 


2.54 


2.56 


2.58 


2.60 


1.76 


1.78 


1.80 


1.82 




B.E.(MeV) 


1.54 


6.59 


13.67 


22.29 


0.92 


4.74 


9.98 


16.34 




M(MeV) 


4119.38 


4114.33 


4107.25 


4098.63 


10739.60 


10735.82 


10730.52 


10724.26 


2 + 


r»(fm) 


2.28 


1.07 


0.75 


0.60 


1.70 


0.75 


0.54 


0.44 




Pi(%) 


99.98 


99.97 


99.97 


99.97 


99.76 


99.79 


99.85 


99.85 




P 2 (%) 


0.00 


0.00 


0.00 


0.00 


0.04 


0.03 


0.02 


0.02 




P 3 (%) 


0.02 


0.03 


0.03 


0.03 


0.20 


0.17 


0.13 


0.10 



to the case of X(3872) l44Tl . With the present approach we 
suggest the state (D^B^) S [J P - + ] to be a candidate of the 
molecule but not a good one due to the strong dependence of 
the binding solutions on the cutoff parameter. 

It is necessary to mention that there exist fourteen chan- 
nels with the quantum numbers S,I(J P ) = 1, A(l + ). It is 
very hard to solve a 14 x 14 matrix Shrodinger equation. 
Due to the large mass gap between the threshold of D*B* s (oi 
D*B*) and that of DB* and the strong repulsive interaction 
coming from the D-wave, we expect the channels Z)*B*( 3 Di), 
D*B* s ( 5 Di), D*B*CDi) and D* S B*{ 5 D X ) to provide negligible 
contributions, which can be clearly seen from the previous 
calculation. Therefore, we omit these four channels in our cal- 
culation, which should be a very good approximation. For the 
state (D^* ) B^* } ) S [J P = 1 + ], it is very interesting that we obtain 
a loosely bound state with a reasonable cutoff 1.24 ~ 1.30 
GeV. If the cutoff is set to be 1.24 GeV, the binding energy 
relative to the threshold of DB* is 2.50 MeV and the root- 
mean-square radius is 1 .45 fm, see Table IVIIII The channel 
DB*( 3 S i) provides the largest contribution, with a probability 
of 58.10%. The total contribution of the D-wave channels is 
very small, less than 0.1%, which also proves that our previ- 
ous approximation is reasonable. With small binding energy 
and large root-mean-square radius and reasonable cutoff, it is 



very natural to interpret the state (D w B w ) S [J P = 1 + ] as a 
molecule. 

For the state (D i * ) B i * ) ) s [J P = 2 + ], if the cutoff parameter is 
fixed to be 2.20 GeV, the binding energy relative to threshold 
of D*B* is 1.23 MeV and the root-mean-square radius is 2.25 
fm. The dominant channel is D*B*( 5 52), with a probability 
of 88.81%. The other S-wave channel D*B*( 5 S 2 ) provides the 
second largest contribution, 10.90%, and the total contribution 
of the D-wave is 0.3%. With such a large cutoff, whether the 
state (D ( * ) B^* ) ) S [J P = 2 + ] exists as a molecule or not is beyond 
the present meson exchange model. Further detailed study 
with other approaches will be helpful to settle this issue. 



C. The Results for The Systems with Strangeness 5=2 

For the systems with strangeness 5=2, there dose not 
exist the long-range pion exchange, but there are mediate- 
range sigma and eta exchanges and the short-range phi ex- 
change. We summarize the numerical results for the systems 
(£>(*)£>(*)) sj and (B ( * ) B ( * ) ) SS in Table US and for the system 
(£)Mfl«) M in TableH 
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TABLE VIII: The numerical results of the D M B ( ' 


system with strangenes 


is S = 1. The binding 


energy of the state (D M B M ) S [J P 


= + ] is 


relative to the DB S threshold while those of the states (D i ' ) B { ' ) ) s [J P = 1 + ] 


and I'D'** R ( * 


tsl J — 


2^ 


] correspond to the thresholds of DB\ and 


D*B* respectively. 














J 






(£><* 




'h 




A(GeV) 


1.28 


1.30 






1.32 


1.34 


B.E.(MeV) 


6.72 


22.10 






43.11 


68.73 


A If f~K If ~\ T\ 

M(MeV) 


7226.81 


7211.43 






7190.42 


7164.80 


r ms (fm) 


0.92 


0.55 






0.43 


0.36 


Pi(%) 


50.10 


36.04 






29.33 


25.30 


P 2 (%) 


25.66 


26.98 






25.07 


23.02 


Pi(%) 


12.03 


18.55 






22.96 


26.05 


P 4 (%) 


12.06 


18.30 






22.56 


25.58 


P s (%) 


0.07 


0.05 






0.04 


0.03 


P 6 (%) 


0.08 


0.06 






0.04 


0.03 


A(GeV) 


1.24 


1.26 






1.28 


1.30 


H 1 — ' /~\ If ~\ J \ 

B.E.(MeV) 


2.50 


14.97 






32.88 


54.86 


H If /~K If ~\ J\ 

M(MeV) 


7280.13 


7267.66 






7249.75 


1221 m 


r rml (fm) 


1.45 


0.63 






0.47 


0.39 


Pi(%) 


19.11 


23.61 






22.14 


20.40 


1+ Pz(%) 
P 3 (%) 


58.10 


35.98 






28.14 


23.97 


4.65 


8.04 






9.71 


10.69 


P 4 (%) 


4.63 


8.16 






9.86 


10.80 


P 5 (%) 


6.73 


11.95 






14.87 


16.84 


Pi(%) 


6.72 


12.21 






15.26 


17.28 


P 7 (%) 


0.00 


0.00 






0.00 


0.00 


Pg(%) 


0.00 


0.00 






0.00 


0.00 


P 9 {%) 


0.04 


0.03 






0.02 


0.01 


Pio(%) 


0.03 


0.02 






0.02 


0.01 


A(GeV) 


2.20 


2.22 






2.24 


2.26 


H 1 — ' /~\ If ~\ T \ 

B.E.(MeV) 


1.23 


5.68 






12.65 


21.29 


H If /~K If ~\ T \ 

M(MeV) 


7422.79 


7418.34 






7411.37 


7402.73 




2.25 


0.95 






0.62 


0.49 


2+ Pd%) 


10.90 


26.64 






36.08 


40.94 


Pi(%) 


88.81 


73.18 






63.83 


59.01 


P 3 (%) 


0.02 


0.01 






0.01 


0.00 


Pa(%) 


0.01 


0.00 






0.00 


0.00 


Ps(%) 


0.18 


0.12 






0.06 


0.03 


P 6 (%) 


0.09 


0.05 






0.02 


0.01 



1. (» W Z) W ).„ and (B W S W ) JS 

If we fix the cutoff parameter between 2.76 GeV and 2.82 
GeV, we obtain a bound state of (D ( *'D ( * ) ). SJ [7 P = + ] with 
binding energy 2.43 ~ 28.53 MeV and root-mean-square ra- 
dius 1.92 ~ 0.55 fm. The dominant channel is D s D s ( 1 Sq), 
with a probability of 94.69% ~ 87.27. The contribution of the 
D-wave channel D* s D*{ 5 Dq) is very small as expected, less 
than 0.1%. In the bottomed sector, we obtain binding solu- 
tions with a smaller but more reasonable cutoff 1 .90 ~ 1 .96 
GeV. As one can easily read off from Table |IX] when we set 
the cutoff to be 1.90 GeV, the binding energy is 2.27 MeV 
and the root-mean-square is 1.17 fm. The dominant chan- 
nel is BsBs^Sq), with a probability of 90.83%. However, 



the D-wave channel B*B*( 5 Dq) provides a negligible contri- 
bution, 0.29%. Based on the binding solutions, the meson 
exchange approach favors (but not strongly) that the state 
(B { * ) B { * ) )[J P = + ] is a molecule candidate. 

For the state (D ( *'D ( * , ) JJ [7 P = 1 + ], when we tune the cutoff 
parameter between 2.62 GeV and 2.68 GeV, we obtain bind- 
ing energy 0.41 ~ 16.57 MeV and root-mean-square radius 
4.64 ~ 0.71 fm. The dominant channel is [D S D* S ] + ( 3 S i), 
with a probability of 99.98% ~ 99.96%. In the corre- 
sponding bottomed case, we obtain a loosely bound state of 
(BB) SS [J P = 1 + ] with binding energy 0.83 - 11.98 MeV and 
root-mean-square radius 1 .95 ~ 0.54 fm when we set the cut- 
off parameter between 1 .82 GeV and 1.88 GeV. Similar to the 
charmed case, the dominant channel is [B S B*] + ( 3 S i), provid- 
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ing a contribution 99.41% ~ 99.29%. In our present meson 
exchange model, the state (B ( * ) B ( * ) ) SJ [7 P = 1 + ] can be viewed 
as a molecule candidate. 

For the J p = 2 + case, the numerical results very similar to 
those of the J p = 1 + case, seellXl 



2. (Z)«flM)„ 

For the state (D ( * ) B ( * ) ) JJ [7 F = + ], we obtain binding en- 
ergy 0.64 ~ 21.72 MeV and root-mean-square radius 3.11 ~ 
0.52 fm with the cutoff parameter fixed between 2.36 GeV 
and 2.42 GeV. The channel D S B S ( Sq) with a probability of 
96.05% ~ 86.20% is the dominant channel. The probabil- 
ity of the channel D*B*( 5 D ) is very small, about 0.03%, see 
Table H 

For the state (D ( * ) B ( * ) ) SS [J P = 1 + ], when we fix the cut- 
off parameter between 2.34 GeV and 2.40 GeV, we obtain 
binding energy 1.47 ~ 23.41 MeV and corresponding root- 
mean-square radius 2.02 ~ 0.50 fm. The dominant channel 
is D s B*( 3 Si), with a probability of 93.37% ~ 83.13%. How- 
ever, the total contribution of the D-wave is very small, less 
than 0.1%, see Table|X] 

Very similar to the J p = + case, we obtain a bound state of 
(Z)«2?M)„[/ P = 2 + ] with binding energy 2.98 ~ 22.29 MeV 
and root-mean-square radius 1.35 ~ 0.51 fm. The dominant 
channel is D*B*( 5 S 2 ), with a probability of 99.98%. Based on 
the present numerical results, whether these three states exist 
as molecules or not are beyond the present meson exchange 
approach. Further detailed study with other approaches will 
be very helpful to settle this issue. 



IV. CONCLUSION 

In the present paper, we investigate the possible molecu- 
lar states composed of two heavy flavor mesons, including 
DWdW, #(*)#(*) and D^fiW w j m strangeness S = 0, 1 and 
2. In our study, we take into account the S-D mixing which 
plays an important role in the formation of the loosely bound 
deuteron, and particularly, the coupled-channel effect in the 
flavor space. 

In order to make clear the role of the long-range pion ex- 
change in the formation of the loosely bound states, we give 
the numerical results with the one-pion-exchange potential for 
the system with strangeness S = 0, as well as the numerical 
results with the one -boson-exchange potential. 

In our study, we notice that the contribution of the D-wave 
channel with large threshold is almost negligible for the sys- 
tem with a large mass gap among the thresholds of different 
channels. We also notice that when the binding energy is com- 
parable to or even smaller than the mass gap, the effect of 
the mass gap will be magnified by the small binding energy, 
which is similar to the X(3872) case Hi . 

In the sector with strangeness 5=0, our results 
favor that D&D&[I(J P ) = 0(1 + )], B^B^[I(J P ) = 
0(l + ),0(2 + ),l(0 + ), 1(1 + ), 1(2 + )] and D { * ) B ( * ) [I(J P ) 
0(0 + ),0(l + )] are good molecule candidates. For these 



states, the long-range pion exchange is strong enough to 
from the loosely bound states, and the mediate-range eta 
and sigma exchanges and the short-range rho and omega 
exchanges are helpful to strengthen the binding. How- 
ever, the possibilities of the states D W D W [/(7 P ) = 0(2 + )] 
and D^*'B^*'[I(J P ) = 1(1 + )] as molecules are ruled out 
by the present meson exchange approach. Whether 
the states D^D < -*\l(J p ) = 1(0 + ), 1(1 + ), l(2 + )]and 
D^B^[I(J P ) = 0(2 + ), 1(0 + ), 1(2 + )] exist as molecules 
or not are beyond the present meson exchange model. Further 
detailed study with other approaches will be very helpful to 
settle this issue. 

In the sector with strangeness S = 1, from our results 
the states_ (D W D W ) ! [/(7 P ) = 1 + ], (B M B ( * ) ) S [J P = 1 + ] 
and (D (t) S w ) S [J P = + , 1 + ] are good candidates of molec- 
ular states. The states (B i * > B { * ) ) s [J p = + ,2 + ] may also 
be molecules, but not strongly favored. Whether the states 
(flWflWU/ = + ,2 + ] and (D M B^) S [J P = 2 + ] exist as 
molecules or not are beyond the present meson exchange 
model, and further study are needed. 

For the system with strangeness 5=2, our results 
favor (but not strongly) that the states (B i * ) B { * ) ) ss [J p = 
+ , l + ,2 + ] exist as molecules. However, with the present 
meson exchange model we cannot draw a definite conclu- 
sion about whether the states (£> ( * ) D ( * ) ). !S |7 /> = + , 1 + , 2 + ] and 
(D i * ) B { * ) ) ss [J p = + , 1 + , 2 + ] are molecules or not, and the fur- 
ther detailed study with other approaches are needed. 

It is very interesting to search the predicted exotic hadronic 
molecular states experimentally. These molecular candidates 
cannot directly fall apart into the corresponding components 
due to the absence of the phase space. For these molecular 
states with double charm, they cannot decay into a double 
charm baryon plus a light baryon. The masses of the light- 
est doubly-charmed baryon and light baryon are 3518 MeV 
and 938 MeV, respectively, corresponding to H* c and proton as 
listed in PDG 14111 . The mass of the molecular state is around 
3850 MeV and much smaller than the sum of the masses of a 
doubly-charmed baryon and a light baryon. Therefore such a 
decay is kinematically forbidden. 

However, the heavy vector meson within the exotic molec- 
ular state may decay. The D* mainly decays to Dn via strong 
interaction. It also decays into Dy. Similarly, £)*, B* and B* 
dominantly decay to D s y, By and B s y via electromagnetic in- 
teraction, respectively. For example, the main decay modes of 
the exotic double-charm molecular state with one D* meson 
are DDy and DDn. If the molecular candidates contain two 
heavy pseudoscalar mesons only, they are stable once pro- 
duced. The D* meson may also decay via weak interaction. 
For these exotic molecules with double bottom or both one 
charm and one bottom, their decay behavior is similar to that 
of the molecular state with double charm. 

The above typical decay modes provide important informa- 
tion to further experimental search. Although very difficult, 
it is still possible to produce such heavy systems with double 
bottom or both one charm and one bottom at LHC. 
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TABLE IX: The numerical results for the Z) < * ) £) ( * ) /S ( * , -B ( * , systems with strangeness S = 2. The binding energy of the state 
(D M D ( * ) ) s J(B i * ) B i ' ) ) ss [J P = + ] corresponds to the D S DJB S B S threshold while that of the state (£» ( * , £> ( * , ) IJ /(5 ( * ) 5 ( * ) )., I corresponds to the 
threshold of D S D'JB S B' S . 



J p 






(£>'* 








(S<« 








A(GeV) 


2.76 


2.78 


2.80 


2.82 


1.90 


1.92 


1.94 


1.96 




B.E.(MeV) 


2.43 


8.68 


17.56 


28.53 


2.27 


6.97 


13.30 


20.99 


r\+ 




M(MeV) 


3934.55 


3928.30 


3919.42 


3908.45 


10730.33 


10725.63 


10719.30 


10711.61 


r r ,„Afm) 


1.92 


1.00 


0.70 


0.55 


1.17 


0.67 


0.49 


0.40 




Pi{%) 


94.69 


91.37 


89.04 


87.27 


90.83 


87.05 


84.67 


83.00 




Pi(%) 


5.29 


8.59 


10.91 


12.68 


8.88 


12.65 


15.05 


16.76 




Pi(%) 


0.02 


0.03 


0.04 


0.05 


0.29 


0.30 


0.28 


0.24 




A(GeV) 


2.62 


2.64 


2.66 


2.68 


1.82 


1.84 


1.86 


1.88 




B.E. (MeV) 


0.41 


3.75 


9.32 


16.57 


0.83 


3.39 


7.16 


11.98 


1 


M (MeV) 


4080.38 


4077.04 


4071.47 


4064.22 


10778.27 


10775.71 


10771.94 


10767.12 


r ms (fm) 


4.64 


1.49 


0.94 


0.71 


1.95 


0.96 


0.68 


0.54 




Pi{%) 


99.98 


99.96 


99.96 


99.96 


99.41 


99.29 


99.33 


99.41 




P 2 (%) 


0.01 


0.02 


0.02 


0.01 


0.24 


0.27 


0.24 


0.20 




P 3 (%) 


0.01 


0.02 


0.02 


0.03 


0.35 


0.44 


0.43 


0.39 




A(GeV) 


2.60 


2.62 


2.64 


2.66 


1.82 


1.84 


1.86 


1.88 




B.E.(MeV) 


0.86 


4.79 


10.81 


18.46 


0.22 


2.48 


6.24 


11.15 




M(MeV) 


4223.74 


4219.81 


4213.79 


4206.14 


10830.61 


10828.32 


10824.62 


10819.72 


2 + 


r rmj (fm) 


3.13 


1.29 


0.86 


0.66 


3.86 


1.10 


0.71 


0.55 




P,(%) 


99.98 


99.97 


99.98 


99.98 


99.75 


99.61 


99.64 


99.70 




Pi(%) 


0.00 


0.00 


0.00 


0.00 


0.04 


0.06 


0.05 


0.04 




P 3 (%) 


0.02 


0.02 


0.02 


0.02 


0.21 


0.33 


0.30 


0.26 



V. ACKNOWLEDGEMENT 

This project was supported by the National Natural Sci- 
ence Foundation of China under Grants 11075004, 11021092, 
11261130311, 11222547, 11175073, 11035006, Ministry of 
Science and Technology of China (2009CB825200), the Min- 
istry of Education of China (FANEDD under Grant No. 



200924, DPFIHE under Grant No. 2009021 1 120029, NCET, 
the Fundamental Research Funds for the Central Universities) 
and the Fok Ying-Tong Education Foundation (No. 131006). 
This work is also supported in part by the DFG and the NSFC 
through funds provided to the sinogermen CRC 110 "Symme- 
tries and the Emergence of Structure in QCD". 



[1] N. A. Tornqvist, (2003), arXiv:hep-ph/0308277 [hep-ph]. 

[2] F. E. Close and P. R. Page, Phys. Lett. B578, 1 19 (2004), 

arXiv:hep-ph/0309253. 
[3] E. S. Swanson, Phys. Lett. B588, 189 (2004), arXiv:hep- 

ph/03 11229. 

[4] N. A. Tornqvist, Phys. Lett. B590, 209 (2004), arXiv:hep- 
ph/0402237. 

[5] X. Liu, X.-Q. Zeng, and X.-Q. Li, Phys. Rev. D72, 054023 

(2005), arXiv:hep-ph/0507177. 
[6] S.-L. Zhu, Phys. Lett. B62S, 212 (2005), arXiv:hep- 

ph/0507025. 

[7] M. T. AlFiky, F. Gabbiani, and A. A. Petrov, Phys. Lett. B640, 

238 (2006), arXiv:hep-ph/0506141 [hep-ph]. 
[8] Y.-R. Liu, X. Liu, W.-Z. Deng, and S.-L. Zhu, Eur. Phys. J. 

C56, 63 (2008), arXiv:0801.3540 [hep-ph]. 
[9] X. Liu, Y.-R. Liu, W.-Z. Deng, and S.-L. Zhu, Phys. Rev. D77, 

034003 (2008), arXiv: 07 11.0494 [hep-ph], 
[10] X. Liu, Y.-R. Liu, W-Z. Deng, and S.-L. Zhu, Phys. Rev. D77, 

094015 (2008), arXiv:0803.1295 [hep-ph]. 



[1 1] C. E. Thomas and F. E. Close, Phys. Rev. D78, 034007 (2008), 

arXiv:0805.3653 [hep-ph]. 
[12] S. H. Lee, K. Morita, and M. Nielsen, Phys. Rev. D78, 076001 

(2008) , arXiv:0808.3168 [hep-ph]. 

[13] F. Close and C. Downum, Phys. Rev. Lett. 102, 242003 (2009), 

arXiv:0905.2687 [hep-ph]. 
[14] G.-J. Ding, Phys. Rev. D79, 014001 (2009), arXiv:0809.4818 

[hep-ph]. 

[15] G.-J. Ding, J.-F. Liu, and M.-L. Yan, Phys. Rev. D79, 054005 

(2009) , arXiv:0901.0426 [hep-ph]. 

[16] X. Liu, Z.-G. Luo, Y.-R. Liu, and S.-L. Zhu, Eur. Phys. J. C61, 

411 (2009), arXiv:0808.0073 [hep-ph]. 
[17] X. Liu and S.-L. Zhu, Phys. Rev. D80, 017502 (2009), 

arXiv:0903.2529 [hep-ph]. 
[18] Y.-R. Liu and Z.-Y. Zhang, Phys. Rev. C79, 035206 (2009). 
[19] S. H. Lee, K. Morita, and M. Nielsen, Nucl. Phys. A815, 29 

(2009), arXiv:0808.0690 [hep-ph]. 
[20] P. Ortega, J. Segovia, D. Entem, and F. Fernandez, AIP Conf. 

Proc. 1257, 331 (2010). 



15 



TABLE X: The numerical results of the D M B (,) system with strangeness 5=2. The binding energy of the state (£> ( *'B ( * ) ) (! [y /> = + ] is relative 
to the threshold of D S B S while that of the state (D ( * ) S ( * ) ) SI [7 P = 1 + ] corresponds to the D S B' S threshold. 



J P (DB) SS 





A(GeV) 


2.36 


2.38 


2.40 


2.42 




B.E.(MeV) 


0.64 


5.16 


12.43 


21.72 


+ 


M(MeV) 


7334.15 


7329.63 


7322.36 


7313.07 


r,.,„ r (fm) 


3.11 


1.06 


0.68 


0.52 




Pi (%) 


96.05 


91.15 


88.25 


86.20 




Po(%") 


3.94 


8.82 


11.72 


13.17 




p 3 {%) 


0.02 


0.03 


0.03 


0.03 




A(GeV) 


2.34 


2.36 


2.38 


2.40 




B.E. (MeV) 


1.47 


6.63 


14.13 


23.41 




M (MeV) 


7382.42 


7377.26 


7369.76 


7360 48 




r^„, r ffm) 


2.02 


0.92 


0.63 


0.50 




Pl(%) 


93.37 


88.46 


85.36 


83.13 




p 2 (%) 


0.00 


0.00 


0.00 


0.00 




p 3 (%) 


3.85 


6.95 


9.07 


10.69 




p 4 (%) 


0.01 


0.01 


0.01 


0.01 




p*(%) 


2.75 


4.55 


5.54 


6.14 






0.00 


0.00 


0.00 


0.00 




p 7 (%) 


0.02 


0.02 


0.02 


0.02 




A(GeV) 


2.26 


2.28 


2.30 


2.32 




B.E. (MeV) 


2.98 


7.96 


14.47 


22.29 




M(MeV) 


7524.72 


7519.74 


7513.23 


7505.41 


2 + 


'•mi (fm) 


1.35 


0.82 


0.62 


0.51 




Pl(%) 


99.93 


99.94 


99.95 


99.96 




P 2 (%) 


0.01 


0.01 


0.01 


0.01 




P 3 (%) 


0.06 


0.05 


0.04 


0.03 



[21] S. Ohkoda, Y. Yamaguchi, S. Yasui, K. Sudoh, and A. Hosaka, 
Phys. Rev. D86, 014004 (2012), arXiv: 111 1.2921 [hep-ph]. 

[22] S. Ohkoda, Y. Yamaguchi, S. Yasui, K. Sudoh, and A. Hosaka, 
(2012), arXiv: 1209.0144 [hep-ph]. 

[23] F. Aceti, R. Molina, and E. Oset, PoS QNP2012, 072 (2012). 

[24] M. B. Voloshin and L. B. Okun, JETP Lett. 23, 333 (1976). 

[25] A. De Rujula, H. Georgi, and S. L. Glashow, Phys. Rev. Lett. 
38, 317 (1977). 

[26] N. A. Tornqvist, Nuovo Cim. A107, 2471 (1994), arXiv:hep- 

ph/93 10225 [hep-ph]. 
[27] N. A. Tornqvist, Z. Phys. C61, 525 (1994), arXiv:hep- 

ph/93 10247. 

[28] M. Mattson et al., SELEX Collaboration, Phys. Rev. Lett. 89, 

112001 (2002), arXiv:hep-ex/0208014 [hep-ex]. 
[29] A. Ocherashvili et al., SELEX Collaboration, Phys. Lett. B628, 

18 (2005), arXiv:hep-ex/0406033 [hep-ex]. 
[30] B. Aubert et al., BABAR Collaboration, Phys. Rev. D74, 

011103 (2006), arXiv:hep-ex/0605075 [hep-ex]. 
[31] R. Chistov et al, BELLE Collaboration, Phys. Rev. Lett. 97, 

162001 (2006), arXiv:hep-ex/0606051 [hep-ex]. 
[32] M. B. Wise, Phys. Rev. D45, R2188 (1992). 
[33] A. F. Falk and M. E. Luke, Phys. Lett. B292, 119 (1992), 

arXiv:hep-ph/9206241 [hep-ph]. 
[34] B. Grinstein, E. E. Jenkins, A. V. Manohar, M. J. Savage, 

and M. B. Wise, Nucl. Phys. B380, 369 (1992), arXiv:hep- 

ph/9204207 [hep-ph]. 
[35] R. Casalbuoni et al., Phys. Rept. 281, 145 (1997), arXiv:hep- 



ph/9605342 [hep-ph]. 
[36] Y.-B. Dai and S.-L. Zhu, Eur. Phys. J. C6, 307 (1999), 

arXiv:hep-ph/9802227 [hep-ph]. 
[37] F. Navarra, M. Nielsen, M. Bracco, M. Chiapparini, and 

C. Schat, Phys. Lett. B489, 319 (2000), arXiv:hep-ph/0005026 

[hep-ph]. 

[38] S. Ahmed et al, CLEO Collaboration, Phys. Rev. Lett. 87, 
251801 (2001), arXiv:hep-ex/0108013 [hep-ex]. 

[39] C. Isola, M. Ladisa, G. Nardulli, and P. Santorelli, Phys. Rev. 
D68, 114001 (2003), arXiv:hep-ph/0307367 [hep-ph]. 

[40] M. Bando, T. Kugo, and K. Yamawaki, Phys. Rept. 164, 217 
(1988). 

[41] K. Nakamura and P. D. Group, Journal of Physics G: Nuclear 

and Particle Physics 37, 075021 (2010). 
[42] I. W. Lee, A. Faessler, T. Gutsche, and V. E. Lyubovitskij, Phys. 

Rev. D80, 094005 (2009), arXiv:0910.1009 [hep-ph]. 
[43] R. Machleidt, K. Holinde, and C. Elster, Phys. Rept. 149. 1, 

(1987). 

[44] N. Li and S.-L. Zhu, Phys. Rev. D86, 074022 (2012), 

arXiv:1207.3954 [hep-ph]. 
[45] W. Meguro, Y.-R. Liu, and M. Oka, Phys. Lett. B704, 547 

(2011), arXiv: 1105.3693 [hep-ph], 
[46] N. Li and S.-L. Zhu, Phys. Rev. D86, 014020 (2012), 

arXiv: 1204.3364 [hep-ph]. 



16 



VI. APPENDIX 

The functions Hi etc are defined as, 



and 



Hq(A, qo, m, r) 
Hi(A,q ,m, r) 
H 3 (A,q ,m, r) 
Mi (A, qo,m,r) 
Mi(A,qo,m,r) 

where, 



u 

4^ 
ir_ 

An 

Yin 

.3 



u Lu 



Y(ur)-^Y(xr) 
u 



r/3 2 
it 



Y(xr) 



2u 3 

t 
2m 3 

-— jl [cos(0r) - e~ xr ] + ^e~ xr 
4n I Br L J 26> 3 



Z{ur) - ^Z(xr) - 



U 



cos (#r) 



3 sin (8r) 3 cos (8r) 



Or 



2 r 2 



1 y 3 yB 2 



j8 2 = A 2 - m 2 , 



2 2 2 

u = m — q , 



e 2 = -(m 2 - AT 2 = A 2 - ^ 



F(jc) = — , Z(jc) = |l + - + 4-| ^M, 

X \ X X L 1 



■\- + \\Y(x), 



Fourier transformation formulas read: 
Ho(A, q () , m, r), 
-Hi(A,q ,m,r), 



u 2 + q 2 



u 2 + q 2 

mi 
u 2 + q 2 



-i/ 3 (A, q , m, r)kjj - ^ #i( A > <7o, m, r)6 u , 



(37) 



where, kjj = 3^ - Sij. 

We summarize the isospin-dependent coefficients in Ta- 
bles IXIIXIllXlTTl and the time component of the transferred 
momentum used in our calculation in Table IXIVl 
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TABLE XI: The isospin-dependent coefficients for the Z) ( * , £) ( * ) and D ( * ) B i * ) systems with strangeness 5 


= 0. The coefficients of the B W B W 


systems, which are not shown, are similar to those of the D^D^*\ 






1 = 


DD TDD*1 


[DD'} + 


D'D' 


1)1) 


3 _a , 1 . .a . _fl 
+ ^ CO + cr 




-\if + |t? c - |p c + \o) c 


yuu j_ 


3 ^d _■_ 1 . .d , ~d 




-^(-!^ s + l 6 i s - |p s + V) 


[DD*] + 




3 „rf , 1 rf 1 „rf 

— 2P + 9 W + " 
+ (-§«*+ i?f - |p 6 + 


i ( -^/ + y-|p/ + X) 
+ 7f(-|^ s + W - ¥ + 1^ 


D*D* 






-\n h + -\p h + \d' +a h 


I - i 


U U [ I s U J_ 


[DD'] + 


D'D* 


1)1) 


1 „(7 , I , ,(2 1 ~-Ct 

-p + -6J + u 




\if + \rf + \p c + i<o° 


rnnn 

yuu j_ 


—(l\ 7r e _i_ I _l I n <? _i_ L,,f \ 




^ + y + i P / + X) 


yuu j + 




1 d , 1 , d , d 

^p + ^co + cr 

./I p I p 1 p 1 p \ 

+(K + + + ) 


+T2 (i 2 ns + ^ s + y* + ^ 


u L) 






\n h + \jf + \f) + \(x) h + (T h 


1 — \) 


UD UD 




D*B* 


n r 

UD 


~jP + 2 




-|^ + |77 c -fp c + 3^ 


DB* 


_3 n rf 4- 4- rr"' 


3 c , I e 3 f , 1 e 

~W + vf - iff + i<s 


-Inf + lrjf-y + W 


ix 5 
I J D 




— 2-p + ^ai + tr 


-\7? + \rf-\f? + \«? 


U D 






-\n h + \rj> - \p h + \d' + a h 


/ = 1 




D*B 


D*B* 


DB 


|p fl + ±w" + cr" 




\lf + \t] c + \p c + \cf 


DB* 


|p rf + iti/* + cr 1 * 


iy + Ijf + \p e + \uf 


W + y + y + x 


D*B 




\p d + {u> d + cr'' 


i«« + + ^p^ + w 


D*B* 






\n h + \rf + ip* + ia** + cr* 



TABLE XII: The isospin-dependent coefficients for the (D {,) D M ) S and (D M B M ) S systems with stranj 
(j} ( *'j}0O) s systems, which are not shown, are similar to those of the (D'* ) £) ( * , ) I . 


;eness 5 = 1. The coefficients of the 




DD S 


DD' S 


D*D S 




D'D' S 


DD, 


cr a + K' a 










+ K c + K' c 


dd: 




<r d + K e + K* e 


-\Tf+K** 


-w 


+ K f + K'f 


D'D S 






a d +K e + K' e 




+ K s + K' s 


d-d; 








-\rj h + a h + K h + K' h 




DB S D S B 


DB] D*B 


D*B S D S B* 


D'B' S 


D*B* 


DB S 


cr" K ,a 








-w 


K c + K' c 


D S B 


<r" 








K c + K' c 


-w 


db: 




cr rf K e + K*" 


-\rf e K* d 


-w 


Kf + K* f 


D*B 




a" 


K"' -\jf 




-w 


D'B S 






cr rf K" + K' e 






D S B* 






<T d 




-w 


d-b; 








+<T h 


K h + K ,h 


D*B* 
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TABLE XIII: The isospin-dependent coefficients for the (D M D < -* ) ) SS and (D M B M ) SS systems with strangeness 5 = 1. The coefficients of the 
(gMg») j( systems, which are not shown, are similar to those of the (D'-* > D'-* ) ) !1S . 

1=1 DsDs [DsD'A- [£>.,£>;]+ DjD] 

D S D S (/>" + a a \rf + 4> c 

[D S D* S ]- <p" + - (y + cf>n x J- { y + J- { y + $ S) 

[D t D*] + ^ + ^ + (|rf + ^) ^(f n f + ¥) + ^(l I s + <P S ) 

D' S D* y- + </>'' + a h 

D S B S D S B* S D* S B S D* S B* S 

D S B S (p a + a a \if + fi 

£>,6* <j> d + <r d \rf + <t> e \rf + $ 

D*B S (f + a d y + cf> s 

D;b: y + </>"+ a" 



TABLE XIV: The time component of the transferred momentum, q , used in our calculation. The other values which are not given are zero. 





DD 






66 








DB 




Process 




9o 


Process 




qo 


Process 


9o 


Process 


qo 



DD S -> D S D 
D S D D* S D* 
D S D D*D* S 
D S D* DD* 
DD* —> D*D 
D S D* -> D*D 
D S D* D*D S 
DD* — > D*D 

DD* S -> D '°s 
DD* D*D* 
I), I) -> I) I) 
I), I) > I) I), 

DD* -> D*D* 
DD* -> D*D* 

DID* d*d: 



m Ds - m D 

(m 2 D , +m 2 D )-(m 2 [)! +m 2 D ,) 



2(m D » +m D * ) 



2(m D « +m D <) 
(m 2 D ,+m 2 Ds )~{m 2 D +m 2 D ,) 
2(m D +m D *) 

m D , - m D 

(m 2 D , +m 2 D , )-(m 2 D ^+m 2 D ) 
2(m D ^+m D *) 

m D , - m Ds 



"Id; - m D 

(m 2 D , +m 2 D , )-(n? Di -m 2 D ) 

2(in Ds +m D *) 
2 2 
m D-- m D 



2(m » +m c ») 

m D> +m \ s - 2 '"D' 

2(m B « +m D ») 
2m 2 D ,-m 2 D -m 2 D , 
2(m D * +m D *) 

m Q.-"'o 
2(m » +m D ») 

m D - - m D > 



BB S — > 6,6 
6,6 -> 5*6* 
6,6 5*5* 
B S B* -> 55* 

BB* -> 5*5 
6,6* -> 5*5 
B S B — > B B s 

BB; -> b* s b 

BB* S -> B*B S 
BB* -> 5*5* 
6,6* -> 5*5* 
6,6* -> 5*5* 
66* -> 5*5* 
66* 5*5* 
5t5* 6*5! 



m Bs - m B 

(m 2 , +m|)-(m| a +m 2 , ) 

2(m B * +/» B *) 
(m^» +m 2 ! )-(m 2 , +m 2 ) 

2(m B > +m B «) 
(m^, +m 2 ! )-(m- g +m 2 ,) 

2(m B +m B *) 

m B > - m B 

(m 2 , +m 2 , )-(m|+m| ) 
2(m«, +»ifi» ) 

m B , - m Bs 



m B > - m B 

(m 2 , +m 2 , )-(m|+m| i ) 
m 2 * 



2(m fl » +m s ») 
m 2 B ,+m 2 Br 2i,r B , 
2(m B - +m s «) 

2(m fi . +ra a «) 
2 2 

2(»i s . +m B . ) 

ot s ? - m B ' 



D S B DB S 
D S B -> D*B* 
DB S D*B* 
D S B* DB] 

DB* D*B 
D S B* 0*6 
0,5* -> 0*5, 

D6J ^ 0*6 
DB* S D*B S 
DB* D*B* 
0,6* ^ £»;6* 
£),6* ^ 0*6* 
06J ^ Z)*6* 
DB* S -> D*B* S 
DIB* d*b; 



2(mD+m Bs ) 
(m^, +m|)-(m^ , +m 2 Ds ) 

2(m £1 » +m B «) 
(m^» +m 2 s )-(m 2 , +m 2 D ) 

2(tn D * +m e *) 
( m D, +m B « H»j«J, ) 

2(m a ,+m B «) 
(m^»-m|,)-(m|,+m|) 

2(m «+m B ) 
("J,, +"i B »)-("iJ, i +m|) 

2(m D » +m B ) 
(m 2 , +m 2 D , )-(m Bs +mj, s ) 
2(m c « +m Bs ) 

2(m c . +m B ) 
(m 2 , +m 2 r )-(m 2 D +m 2 ^ 

2(m +m B .) 
2 2 S 

'V-'"b 

2(m D « +;« B . ) 

"'Dr"'o* 
2(m D » +m B «) 

2(m£,« +m B «) 
(m^, +m^» )-(mJ+m^» ) 
2(m D « +m B ») 

2 2 



2(m£,«+m B .) 
(<. +m O.)-( m B . +m fl«) 



2(/71d* +m B * ) 



DB S 
D S B - 
DB S - 
D* S B- 
D*B 

d;b- 

D* S B 
D*B S 
D*B S 
D*B- 
D* S B- 
D*B- 
D*B S 
D*B S 
D*B* 



. D S B 
D*B* S 
D*B* S 
D*B S 

> DB* 
D S B* 

> DB* S 

> D S B* 
>DB' S 

D*B* 
D* S B* 
D*B* S 
>D* S B* 
>D*B* S 

> D*B* 



( m B , +m b,'~''"i! + " ! Z>> 

2(M[)+m B ) 
{m 2 B , +m 2 D ^)-(m- D , +m|) 

2(m B » +m B ») 
(m 2 D , +m| s )-(m 2 » +mj,) 
2(m B » +m B «) 

2(m . +m Bs ) 
(m 2 D , -m 2 t )-(m 2 D +m 2 ) 

2(m »+m B ) 
(f 2 B ,+m 2 D , )-(m 2 B +m 2 D ) 

2(m D » +;n B ) 

2(m . +;« B ) 
(m^, +m^, s +m 2 D s ) 

2(m ,+m B ») 
(m^, +m 2 , ) 

2(mD+m B *) 
mL —hit. 



2(m £) « +m B » ) 

"'!»-"'! 
2(m D » +m B ») 

(™ D . +m B » )-(»i D « 
2(m e » +m B «) 

(m^, +m| s )-(mJ, +m 2 , ) 
2(m D » +m B ») 

/lit,* -HIT, 



2(m D - +m B 



2(m £) * +/ii B * ) 



